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A systematic study is presented for centrality, transverse momentum (pr) and pseudorapidity (77) 
dependence of the inclusive charged hadron elliptic flow (^2) at midrapidity (I77I < 1.0) in Au+Au 
collisions at ^smn = 7.7, 11.5, 19.6, 27 and 39 GeV. The results obtained with different methods, 
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including correlations with the event plane reconstructed in a region separated by a large pseu- 
dorapidity gap and 4-particle cumulants (w2{4}), are presented in order to investigate non-flow 
correlations and V2 fluctuations. We observe that the difference between V2{2} and V2{4} is smaller 
at the lower collision energies. Values of V2, scaled by the initial coordinate space eccentricity, V2/e, 
as a function of pr are larger in more central collisions, suggesting stronger collective flow develops 
in more central collisions, similar to the results at higher collision energies. These results are com- 
pared to measurements at higher energies at the Relativistic Heavy Ion Collider (\/sjvjv ~ 62.4 and 
200 GeV) and at the Large Hadron Collider (Pb -I- Pb coUisions at = 2.76 TeV). The «2(pt) 

values for fixed pr rise with increasing collision energy within the pr range studied (< 2 CeV/c). A 
comparison to viscous hydrodynamic simulations is made to potentially help understand the energy 
dependence of V2{pt)- We also compare the V2 results to UrQMD and AMPT transport model 
calculations, and physics implications on the dominance of partonic versus hadronic phases in the 
system created at Beam Energy Scan (BES) energies are discussed. 

PACS numbers: 25.75.Ld, 25.75.Dw 



I. INTRODUCTION 

Azimuthal anisotropics of particle distributions rela- 
tive to the reaction plane (plane subtended by the impact 
parameter and beam direction) in high energy heavy- 
ion collisions have been used to characterize the colli- 
sion dynamics [H-Q. In a picture of hydrodynamic ex- 
pansion of the system formed in the collisions, these 
anisotropies are expected to arise due to initial pres- 
sure gradients and subsequent interactions of the con- 
stituents Specifically, differential measurements 
[igj of azimuthal anisotropy have been found to be sensi- 
tive to (a) the equation of state (EOS), (b) thermaliza- 
tion, (c) transport coefficients of the medium, and (d) ini- 
tial conditions in the heavy-ion collisions. Hence it is im- 
portant to study the dependence of azimuthal anisotropy 
as a function of several variables, for example center-of- 
mass energy {-^snn ), collision centrality, transverse mo- 
mentum (pt), and pseudorapidity (77). 

Recently a beam-energy scan (BES) program has be- 
gun at RHIC to study the QCD phase diagram [2^. The 
BES program extends the baryonic chemical potential 
ifiB) reach of RHIC from 20 to about 400 MeV 
The baryon chemical potential decreases with the de- 
crease in the beam energy while the chemical freeze-out 
temperature increases with increase in beam energy p3| . 
This allows one to study azimuthal anisotropy at midra- 
pidity with varying net-baryon densities. Lattice QCD 
calculations suggest that the quark-hadron transition is 
a crossover for high temperature (T) systems with small 
/is or high ^ysNN [13 ■ Several model calculations sug- 
gest that at larger values of fiB or lower s/snn the transi- 
tion is expected to be first order [25l - l27j . Theoretical cal- 
culations suggest a non-monotonic behavior of V2 could 
be observed around this "softest point of the EOS" [28| . 
The softest point of the EOS is usually referred to as 
the temparature/time during which the velocity of sound 
has a minimum value (or reduction in the pressure of 
the system) during the evolution. Non-monotonic vari- 
ation of azimuthal anisotropy as a function of collision 
centrality and -ysNN could indicate the softest point of 
the EOS in heavy- ion reactions [l^. Further it has been 
argued that the observation of saturation of differential 



azimuthal anisotropies V2{pt) of charged hadrons in Au 
+ Au collisions in the s/snn range of 62.4 - 200 GeV is a 
signature of a mixed phase [la] . The new data presented 
in this paper shows to what extent such a saturation ef- 
fect is observed. 

Several analysis methods for V2 have been pro- 
posed [sol - Is^ . These are found to be sensitive in varying 
degrees to non-flow contributions (e.g. correlations due 
to jets, resonances, etc.) and flow fluctuations. V2 mea- 
surements from various methods have been judiciously 
used to constrain these contributions, in addition to pro- 
viding estimates of systematic errors associated with the 
measurements [1^. This is particularly useful for inter- 
preting results of identified hadron V2 values where, due 
to limitations of event statistics, it is not possible to use 
all methods for V2 analysis. The measurements over a 
range of energies may provide insights to the evolution 
of non-flow and flow fluctuations as a function of collision 
energy. 

Inclusive charged hadron elliptic flow measurements at 
top RHIC energies have been one of the most widely 
studied observables from the theoretical perspective. It 
has been shown that transport models, which provide 
a microscopic description of the early and late non- 
equilibrium stages of the system, significantly underpre- 
dict V2 at top RHIC energies, while the inclusion of 
partonic effects provides a more satisfactory explana- 
tion [36| . The new data discussed here will provide an 
opportunity to study the contribution of partonic matter 
and hadronic matter to the V2 measurements as a func- 
tion of ^s/v/v or (T, /^b) by comparisons with models. 

In this paper we present measurements of the second 
harmonic azimuthal anisotropy using data taken in the 
BES program from ~^snn = 7.7 to 39 GeV. We discuss 
the detectors used in the analysis, data selections and 
methods used to determine inclusive charged hadron V2 
in Sections II and HI. Section IV gives V2 results for in- 
clusive charged hadrons from different analysis methods. 
We discuss the centrality, 77, px and ^/snn dependence of 
V2 in Section V, and compare to calculations from trans- 
port models. Finally, a summary of the analysis is pre- 
sented in Section VI. 
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II. EXPERIMENTS AND DATA SETS 

A. STAR detector 

The results presented here are based on data collected 
during the tenth and eleventh RHIC runs (2010 and 
2011) with the STAR detector using minimum-bias trig- 
gers (requiring a combination of signals from the Beam- 
Beam Counters (BBC) [s^l, Zero Degree Calorimeters 
(ZDC) [H, and Vertex Position Detectors (VPD) [H). 
For the 7.7 and 11.5 GeV data, at least one hit in the 
full barrel Time-Of-Flight detector [i^l was required in 
order to further reduce the background. The main Time 
Projection Chamber (TPC) [4l| and two Forward Time 
Projection Chambers (FTPCs) [l^] were used for par- 
ticle tracking in the central region {\ri\ < 1.0) and for- 
ward regions (2.5 < \t]\ < 4.0) respectively. Both the 
TPC and FTPCs provided azimuthal acceptance over 
2tt. The BBC detector subsystem consists of two detec- 
tors mounted around the beam pipe, each located outside 
the STAR magnet pole-tip at opposite ends of the TPC 
approximately 375 cm from the center of the nominal 
interaction point. Each BBC detector consists of hexag- 
onal scintillator tiles arranged in four concentric rings 
that provided full azimuthal coverage. The inner tiles of 
the BBCs, with a pseudorapidity range of 3.8 < jr/l < 5.2 
were used to reconstruct the event plane in one elliptic 
flow analysis. 



B. Event and track selection 

Events for analysis are selected based on collision ver- 
tex positions within 2 cm of the beam axis to reduce 
contributions from beam-gas and beam-pipe (at a radius 
of 4 cm) interactions, and within a limited distance from 
the center of the detector along the beam direction (±70 
cm for the 7.7 GeV data set, ±50 cm for the 11.5 GeV 
data set, and ±40 cm for the 19.6, 27 and 39 GeV data 
sets) . These values are chosen to reduce systematics due 
to variance in detector performance over [yy] < 1.0 while 
retaining sufficient statistics. After quality cuts, about 
4 million - 80% central events remain for 7.7 GeV, 11 
million for 11.5 GeV, 20 mihion for 19.6 GeV, 40 million 
for 27 GeV and 120 miUion for 39 GeV data sets. The 
results from more peripheral collisions are not presented 
due to trigger inefficiencies at low multiplicity. The cen- 
trality was defined using the number of charged tracks 
with quality cuts similar to those in Ref. [l^]. The de- 
tails of the ccntrality determination will be discussed in 
subsection C. The - 80% central events for V2 analysis 
of charged hadrons are divided into nine centrality bins: 
- 5%, 5 - 10%, 10 - 20%, 20 - 30%, 30 - 40%, 40 - 50%, 
50 - 60%, 60 - 70% and 70 - 80%. 

A variety of track quality cuts are used to select good 
charged particle tracks reconstructed using information 
from the TPC or FTPCs. The distance of closest ap- 
proach (DCA) of the track to the primary vertex is taken 



to be less than 2 cm. We require that the TPC and FT- 
PCs have a number of fit points used for reconstruction 
of the tracks to be > 15 and > 5, respectively. For the 
TPC and FTPCs the ratio of the number of fit points to 
maximum possible hits is > 0.52. An additional trans- 
verse momentum cut (0.2 < pT < 2 GeV/c) is applied to 
the charged tracks for the TPC and FTPC event plane 
determination. 



C. Centrality determination 



The centrality classes are defined based on the uncor- 
rected charged particle multiplicity (A^ch*) distribution 
in the TPC for pseudorapidity \ri\ < 0.5 and full azimuth. 



Figure [T] shows the distribution for charged par- 

ticles from the data at ^smn = 7.7, 11.5, 19.6, 27 and 39 
GeV compared to those from Monte Carlo (MC) Glauber 
simulations. The detailed procedures to obtain the simu- 
lated multiplicity are similar to that described in [i^l ■ A 
two-component model [i^l is used to calculate the simu- 
lated multiplicity distribution given by 



dN, 



ch 



drj 



r;=0 



(l-.T)%^ + xA^con 



(1) 



where TVpart is the number of participant nucleons and 
-^coU is the number of binary nucleon-nucleon collisions 
in the simulations. The fitting parameter Upp is the aver- 
age multiplicity per unit of pseudorapidity in minimum- 
bias p + p collisions and x is the fraction of production 
from the hard component. The inelastic nucleon-nucleon 
cross section cr^^fi is extracted from fitting the results 
of available data for total and elastic p + p cross sec- 
tions from the Particle Data Group [431 • The x value is 
fixed at 0.12 ± 0.02 based on the linear interpolation of 
the PHOBOS results at = 19.6 and 200 GeV 

Systematic errors on Upp are evaluated by varying both 
Upp and X within the quoted x uncertainty to determine 
the minimum to describe the data. Since the n^p and 
X are anti-correlated, lower (higher) Upp is used for higher 
(lower) X for systematic error evaluations on A^part- Ta- 
ble U summarizes the parameters in the two-component 
model and cr™^^ in the MC Glauber simulations. The 
cvent-by-cvent multiplicity fiuctuations are included us- 
ing negative binomial distributions [43| . The centrality 
classes arc defined by the fractions of geometrical cross 
section from the simulated multiplicity distributions. For 
each centrality bin, average quantities are calculated in 
the MC Glauber simulations for (Apart), (-/Vcoii), reaction 
plane eccentricity (eRp), participant eccentricity (Cpart), 
root-mean-square participant eccentricity epart{2}, and 
transverse area (S'part)- Eccentricity and transverse area 
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FIG. 1: (Color online) Distribution of uncorrected multiplicity A'^^"^ measured within < 0.5 in the TPC from ^/snn = 7.7 
to 39 GeV in Au + Au collisions shown as black points. The red curves show the multiplicity distributions at ^/sNN ~ 7.7 to 
39 GeV from MC Glauber simulations. See texts for more details about simulations. 



are defined by 



energies. 



^part 
Impart 



(2) 



9 9 ' 



-part 



{2} = yi^IJ- (3) 



xy 



{xy} 



{y'}-{yr 



(4) 
(5) 

(6) 



where the curly brackets denote the average over all par- 
ticipants per event, and x and y are the positions of 
participant nucleons. Systematic uncertainties on those 
quantities arc evaluated by varying parameters for the 
two-component model and by varying the input parame- 
ters in the MC Glauber model. The quoted errors are the 
quadratic sum of the individual systematic uncertainties. 
Table HIl summarizes the centrality classes as well as the 



TABLE I: Summary of Upp and crj^^v with systematic uncer- 
tainties at ^/sjvTv = 7.7, 11.5, 19.6, 27 and 39 GeV. x is set to 
0.12 ± 0.02 for all collision energies. 



y/SNN (GeV) 




o'Tn (mb) 


7.7 


0.89 ± 0.04 


30.8 ± 1.0 


11.5 


1.07 ± 0.05 


31.2 ± 1.0 


19.6 


1.29 ± 0.05 


32.0 ± 1.0 


27 


1.39 ± 0.06 


33.0 ± 1.0 


39 


1.52 ± 0.08 


34.0 ± 1.0 



results obtained by MC Glauber simulations at the five 



III. ELLIPTIC FLOW METHODS 
A. The event plane method 

The event plane method js^l correlates each particle 
with the event plane determined from the full event mi- 
nus the particle of interest, which can be done for each 
harmonic. For any Fourier harmonic, n, the event flow 
vector [Qn) and the event plane angle (\E'„) are defined 
by M 



Qn cos n*,! = Qr, 



QnSVdn'^n = Qny 



^ Wi COS ncf>i 

i 

Y^Wi sinn(/ii, 



-1 Qn 



tan ^ -r-^ j /n, 



(7) 



(8) 



(9) 



where sums extend over all particles i used in the event 
plane calculation, and <f>i and Wi are the laboratory az- 
imuthal angle and the weight for the z*^ particle, respec- 
tively. The reaction plane azimuthal distribution should 
be isotropic or flat in the laboratory frame if the detectors 
have ideal acceptance. Since the detectors usually have 
non-uniform acceptance, a procedure for flatte ning the 
laboratory event plane distribution is necessary |47l l48l|. 

As shown in Eq. pi7|) . the observed is calculated 
with respect to the reconstructed event plane angle '^n 
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where n equals 2 when we use the second harmonic event 
plane and n equals 1 when we use the first harmonic event 
plane. 



,,obs 
^2 



(cos[2(0 -*„)]) 



(10) 



The angular brackets indicate an average over all parti- 
cles in all events. However, tracks used for the V2 calcula- 
tion arc excluded from the calculation of the flow vector 
to remove self-correlation effects. Because the estimated 
reaction plane fluctuates due to finite number of parti- 
cles, one has to correct for this smearing by dividing the 
observed correlation by the event-plane resolution (the 
denominator in Eq. which is the correlation of the 

event plane with the reaction plane. 

^ (cos[2(vI/„ -*,)]) ^^^> 

Since the reaction plane is unknown, the denominator 
in Eq. ([TT|) could not be calculated directly. As shown in 
Eq. p2p . we estimate the event plane resolution by the 
correlation between the azimuthal angles of two subset 
groups of tracks, called sub-events A and B. In Eq. (fT2|) 
C is a factor calculated from the known multiplicity de- 
pendence of the resolution [soj . 

(cos[2(M'„-vI/,,)]) = C^/^MX^I^W)]) (12) 

Random sub-events are used for TPC event plane, while 
pseudorapidity sub-events are used for FTPC/BBC event 
plane. 



1. TPC event plane 

The TPC event plane means the event plane recon- 
structed from tracks recorded by the TPC. For this event 
plane the cj) weight method is an effective way to flatten 
the azimuthal distribution for removing detector accep- 
tance bias. These weights arc generated by inverting the 
(j) distributions of detected tracks for a large event sam- 
ple. The (j) weights are folded into the weight Wi in Eq. ([7]) 
and Eq. ®. 

The re-centering correction [13, is another method 
to calibrate the event plane. In this method, one sub- 
tracts from the Q- vector of each event the Q-vector av- 
eraged over many events. For both the 4> weight and 
re-centering methods, the corrections arc applied in each 
centrality bin, in 2 bins of the primary vertex position 
along the longitudinal beam direction (Vz), and in 2 bins 
for positive/negative pseudorapidity. These corrections 
are determined as a function of data collection time. The 
difference in the effects on V2 from the different flattening 
techniques is negligible. 



2. FTPC event plane 



large acceptance losses from hardware faults caused sig- 
nificant gaps in the azimuthal angle distribution of these 
tracks, preventing use of the weight method because 
of the inability to define </> weights in regions of zero ac- 
ceptance. Thus, only the re-centering method is used for 
the FTPC. 



3. BBC event plane 

In this method the first-order event plane is recon- 
structed using particle trajectories determined from hits 
in the BBC detectors. In this case, 4)i denotes the fixed 
azimuthal angle of the center of the i**^ BBC tile in Eq. ([7]) 
and ([H, and Wi is the fraction of BBC-observed energy 
deposition recorded in tile i: 



Wi 



HA, 



(13) 



The BBC event plane obtained from one BBC detector 
is called a sub-event. A combination of the sub-event 
plane vectors for both BBC detectors provides the full 
event plane. 



V2{BBC} 



(cos[2(0-vl/i)]) 
CV(cos[2(vl/f -vfrf)]) 



(14) 



where C is the constant in Eq. ([T 



1 1 



are sub- 



event plane angles from each BBC detector and ^'i is the 
full event plane angle from both sub-events combined. 

The detector acceptance bias is removed by applying 
the shift method [43. Equation ([T5|) shows the formula 
for the shift correction. The averages in Eq. ([T5|) are 
taken from a large sample of events. In this analysis, the 
correction is done up to the twentieth harmonic. The 
distributions of and ^'f are separately fiattened and 
then the full-event event plane distribution is fiattened. 
Accordingly, the observed and resolution are calcu- 
lated using the shifted (sub)event plane azimuthal angles. 



= * + E„ -[-(sin(2n*))cos(2n*) 



(15) 



More details for the BBC event plane have been described 
in Ref. [H. 



B. The r\ sub-event method 

The 77 sub-event method is similar to the event plane 
method, except one defines the flow vector for each par- 
ticle based on particles measured in the opposite hemi- 
sphere in pseudorapidity: 



Forward-going tracks reconstructed in the two FTPCs 
can also be used to determine the event plane. However, 
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Here ?;2{EtaSubs} denotes the results of the rj sub- 
event method and ^'2,r)+(^2,r;_) is the second harmonic 
event plane angle determined by particles with positive 
(negative) pseudorapidity. An t] gap of \r]\ < 0.075 is used 
between negative (positive) rj sub-event to reduce non- 
flow correlations between the two ensembles. 



C. The cumulant method 

The advantage of the cumulant method is that the 
multi-particle cumulant is a higher-order multi-particle 
correlation formalism which removes the contribution of 
non-flow correlations from lower-order correlations [H, 
Issj . The measured 2-particle correlations can be ex- 
pressed with flow and non-flow components: 

2 e \ ' / 

= V^ + Sn 

Here n is the harmonic number and 6n denotes the non- 
flow contribution. The average should be taken for all 
pairs of particles in a certain rapidity and transverse mo- 
mentum region, and for all events of a data sample. The 
measured 4-particle correlations can be expressed as: 

{e'r.{4>i+,P2-4>3-M^ = vf, + 2-2-vl6n+2Sl (18) 

Thus the flow contribution can be obtained by subtract- 
ing the 2-particle correlation from the 4-particle correla- 
tion: 



((e 



m(0i +02-0,3-04) 



m(0i+02-03-04) 



(19) 



where ((...)) is used for the cumulant. The cumulant of 



order two is just ((e 



m(0i-02) 



»=(' 



,in(0i-02)\ 



1. The cumulant method with generating function 

The GF-cumulant method is computed from a gener- 
ating function [ssj : 



Gn{z) 



M 

n[i 



M 



( * 

\z e 



in(pj 



ze 



')] (20) 



Here z is an arbitrary complex number, z* denotes its 
complex conjugate, M denotes the multiplicity in each 
event, and Wj is the weight (transverse momentum, ra- 
pidity etc.) used in the analysis. The event- wise averaged 
generating function then can be expanded in powers of 
z and z* where the coefficients of expansion yield the 
correlations of interest: 

(G„(z)) = l + 2(e-""^i) + -^*(e"'^'> + 

Il^l(^2li^g-in{<l>i+'p2)^ + £li^g»n(</'l+02)^ (21) 

+.zz*(e™('^i-'^2) )) + ... 

These correlations can be used to construct the cumu- 
lants. More details for the analysis of STAR data have 
been described in Rcf. 11 01. 



2. The Q-cumulants method 

The Q-cumulants method [sof is a recent method 
to calculate cumulants without using nested loops over 
tracks and without generating functions [s^. The ad- 
vantage is that it provides fast (one loop over data) and 
exact non-biased (no approximations and no interference 
between different harmonics) estimates of the correlators 
compared to the generating function cumulants. The cu- 
mulants are expressed in terms of the moments of the 
magnitude of the corresponding flow vector Qn 



M 



(22) 



The single-event average two- and four-particle azimuthal 
correlations can be then formulated as: 



(2) 



M(M-l) 



(4) 



_ \Qn\^+\Q2n\''-2-^[Q2nQ*nQn] 

M{M -1){M -2){M -3) 
^ 2(Af-2).|Qn|'-Ai'(Af-3) 
M{M -1){M -2){M -3) 

The average over all events can be performed as: 

_ i:evcnt(W^(2)),(2). 
Zevont (^(2))j 



(23) 



(24) 



(25) 



■in(0i+02-03-04)^^ 

(W^<4)),(4). 



^events 



(26) 



while the weights are the number of two- and four- 
particle combinations: 

W^2) = M{M-1), (27) 



Wt 



(4) 



M(M-l)(Af-2)(Af-3). (28) 



Choosing the multiplicity weights above can make the 
final mult i- particle azimuthal correlations free of mul- 
tiplicity fluctuations |5l| . However, one can also use 



unit weights treating events with different multiplicity 
equally. The two- and four-particle cumulants without 
detector bias then can be formulated as: 



.{2} 



.{4} 



■2 -((2); 



(29) 



(30) 
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FIG. 2: (Color online) The event plane resolutions for Au + Au collisions at ^/sNN = 7.7, 11.5, 19.6, 27 and 39 GeV as a 
function of collision centrality. Panel (a) shows the resolution of the second harmonic event plane from the TPG < 1). Panel 
(b) shows the resolution for second harmonic event plane from the FTPCs (2.5 < {rjl < 4.0) for 39 GeV and second harmonic 
event plane resolution correction using the first-order event plane from the BBCs (3.8 < < 5.2) for 7.7, 11.5, 19.6 and 27 GeV. 



The reference flow (e.g. integrated over px) can be 
estimated both from two- and four-particle cumulants: 



Vn{2} = V'c„{2} 



(31) 



The event average can be obtained as follows: 

T,cvcnts{w(2'))i{2')i 



in 



E»=l(w(2')) 



(37) 



f„{4} = V-c„{4} 



(32) 



Once the reference flow is estimated, we proceed to the 
calculation of differential flow (e.g. as a function of pr) 
of the particle of interest (POI) which needs another two 
vectors p and q. Particles used to estimate reference flow 
are called reference particles (REP). For particles labeled 
as POI: 



For particles labeled as both POI and REP: 



e 



(33) 



(34) 



Then the reduced single-event average two- and four- 
particle correlations are: 



(2') 



Pn 



nipM - niq 



(4') =[PnQnQ*nQn 
- 2 ■ Mpr^Ql 



nQn PnQnQ2'n 

-2-?«g|0„|2 + 7-g„Q* 



- QnQn + (l2nQ*2n + ' PnQn + 2 ' rUqM 

- 6 • m,]/[(TOpM - 3m,)(M - 1)(M - 2)] 



(35) 



(36) 



,u Eevents(w^(4'>)»(4')» 



Multiplicity weights are: 



U)(2') = VTLpM - niq 



y(4,) E {mpM -mq){M -1){M -2) 



(38) 



(39) 



(40) 



The two- and four-particle differential cumulants without 
detector bias are given by: 



dn{2} = m 



d„{4} = ((4'))-2-((2'))((2)) 



(41) 



(42) 



Equations for the case of detectors without uniform ac- 
ceptance can be found in Ref . [s^l ■ Estimations of differ- 
ential flow are expressed as: 



<{2} 



<{4} 



rfn{2} 
Vcn{2} 

rfn{4} 

-c„{2}3/4 



(43) 



(44) 
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FIG. 3: (Color online) The comparison of V2 as a function of pr between GF-cumulant (open symbols) and Q-cumulant (full 
symbols) methods in Au+Au collisions at ^/snn = 39 GeV. i'2{4} fails in most central (0 - 5%) collisions due to the small 
values of 112 and large V2 flucuations. 
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FIG. 4: (Color online) The px (> 0.2 GeV/c) and rj {\ri\ < 1) integrated V2 as a function of collision centrality for Au + Au 
collisions at = 7.7 GeV (al), 11.5 GeV (bl), 19.6 GeV (cl) , 27 GeV (dl) and 39 GeV (el). The results in the top 

panels are presented for several methods of obtaining «2. The bottom panels show the ratio of V2 obtained using the various 
techniques, with respect to V2{2}. The error bars shown are statistical. 



IV. RESULTS 
A. The event plane resolution 

To investigate the non-flow correlations and V2 fluc- 
tuations of the V2 measurements, the event planes from 



different detectors and the cumulant method are used 
in the analysis. The event planes are determined from 
the TPC in the midrapidity region, and the FTPC/BBC 
at forward rapidity. The rj gap between FTPC/BBC to 
TPC could reduce the non-flow contribution in the V2 
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FIG. 5: (Color online) The V2 as a function of pr for 20 - 30% central An + Au collisions at midrapidity for ^/snn = 7.7 GeV 
(al), 11.5 GeV (bl), 19.6 GeV (cl) , 27 GeV (dl) and 39 GeV (el). The top panels show V2 vs. pr using various methods as 
labeled in the figure and discussed in the text. The bottom panels show the ratio of V2 measured using the various methods 
with respect to 112 {2}. 



measurement [T^ . Figure [2] shows the event plane res- 
olution from TPC (panel (a)) and BBC (FTPC) (panel 
(b)). The resolution of the TPC second harmonic event 
plane increases as the collision energy increases, as the 
resolution depends on the multiplicity and the V2 sig- 
nal [131 • Due to limited statistics, the FTPC event plane 
is used only for the 39 GeV data set where the BBC 
event plane cannot be used because of the poor resolu- 
tion. The resolution of the FTPC event plane is about 
four times lowers than the TPC event plane. The BBC 
is used to determine the event plane for the 7.7, 11.5, 
19.6 and 27 GeV data sets. Note the BBC event plane is 
based on the first harmonic, as the f 1 signal is significant 
in the rapidity region covered by the BBC. The quali- 
tively different centrality dependence of the FTPC and 
BBC event plane resolutions is because of the different 
centrality dependence of vi and W2. 



B. Method comparison 

The comparison of V2 as a function of px between 
the GF-cumulant and Q-cumulant methods is shown in 
Fig. [3] for six collision centralities in AuH-Au collisions 
at ^snn — 39 GeV. The GF-cumulant and Q-cumulant 
methods agree within 5% at all five collision energies. 
Compared to GF-cumulant method, the recently devel- 
oped Q-cumulant is the exact cumulant method [sof . 
The observation of consistency between the two meth- 
ods at BES energies implies the GF-cumulant is a good 



approximation. The cumulant method (GF-cumulant or 
Q-cumulant) used in the analysis does not cause differ- 
ence in the comparison with other experimental results 
and theoretical calculations. To be consistent with the 
previous STAR results, we will hereafter show only re- 
sults from the GF-cumulant method. 

Other method comparisons are shown in Figs. |3] and 
[5] for inclusive charged hadrons in Au -I- Au collisions at 
= 7.7 GeV (al), 11.5 GeV (bl), 19.6 GeV (cl), 
27 GeV (dl) and 39 GeV (el). As the V2 measurements 
from various methods are obtained using charged tracks 
recorded at midrapidity {\ri\ < 1), the statistical errors 
on the results from the different V2 methods are thus 
correlated. The conclusions on the differences in V2 val- 
ues from different methods are based on the systematic 
trends observed for the corresponding ratios with respect 
to V2{2}. Figure m shows V2 integrated over 0.2 < pr < 2.0 
GeV/c and |?7| < 1 versus centrality. For comparison pur- 
poses, the integrated V2 values for all methods are divided 
by the values of the 2-particle cumulant method {v2{2}) 
and plotted in panels (a2) through (e2). The results of 
the 4-particlc cumulants are systematically lower than 
the other methods, except for -!;2{FTPC/BBC}. The 
difference is about 10 - 20% in 39, 27 and 19.6 GeV, 
10 - 15% in 11.5 GeV and 5 - 10% in 7.7 GeV. The 77 
sub-event values for peripheral collisions (50 - 60% to 
70-80%) drop below the 2-particle and TPC event plane 
results, indicating the 77 sub-event method could reduce 
some non-flow correlations for peripheral collisions. Non- 
flow correlations are defined as correlations not related to 
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FIG. 6: (Color online) The 112 over e (Glauber) as a function of pt for various collision centralities (10-20%, 30- 40% and 
50-60%) in Au + Au collisions at midrapidity. Panel (a), (b), (c), (d) and (e) show the results for ^/Inn = 7.7, 11.5, 19.6, 
27 and 39 GeV respectively. The data are from U2{EtaSubs}. The error bars and shaded boxes represent the statistical and 
systematic uncertainties respectively, as described in Sec. IV C. 



the reaction plane. The dominant non-flow correlations 
originating from two-particle correlations (such as HBT 
correlations, resonance decay) scale as 1/N [s^l, where 
N is the multiplicity of particles used to determine the 
event plane. Thus the non-flow contribution is larger 
in peripheral collisions. In mid-central and peripheral 
collisions (10 - 20% to 40 - 50%), the data of 1^2 {BBC} 
from 7.7, 11.5, 19.6 and 27 GeV are consistent with V2{4} 
and lower than other methods. It suggests the first-order 
(BBC) event plane suppresses the second-order non-flow 
and/or fluctuation effects. Within statistical errors, the 
results of W2{FTPC} from Au -I- Au collisions at ^/snn 
= 39 GeV are close to V2{2}, i;2{EP} and W2{EtaSubs} 
in semi-central collisions (10 - 20% to 20 - 30%). In the 
peripheral collisions (30 - 40% to 60 - 70%), z;2{FTPC} 
falls between i;2{EtaSubs} and W2{4}. It indicates that 
the ?7 gap between TPC and FTPC reduces the non-flow 
contribution. 

The pt differential 1)2 from various methods for the 
20 - 30% centrality bin are shown in the upper panels 
of Fig. [SJ For comparison, the V2 from other meth- 
ods are divided by the results of the 2-particle cumu- 
lant method and shown in the lower panels of Fig. [S] 
It can be seen that the difference of W2{2} compared to 
z;2{FTPC/BBC}, V2{2} and z;2{EtaSubs} depends on the 
Pt range. A larger difference can be observed in the 
low Pt region {pT < 1 GcV/c). Beyond = 1 GeV/c 
the difference stays constant in the measured pt range. 
The difference between W2{FTPC/BBC} and t;2{4} is rel- 



atively small and less dependent on pT- It suggests the 
non-flow contribution to the event plane and 2-particle 
correlation methods depends on pT- Based on the inter- 
pretation in Ref. [l|, the difference between V2{2}^ and 
i'2{4}^ is approximately equal to non-flow plus two times 
V2 fluctuations. The fact that the ratio of W2{4} to f2{2} 
is closer to 1 at the lower collision energies indicates the 
non-flow and/or V2 fluctuations in the 1)2 measurement 
depend on the collision energy. One possible explanation 
is that the non-flow correlations from jets presumably 
decrease as the collision energy decreases. The results of 
W2{BBC} are found to be consistent with i'2{4} in 7.7, 
11.5, 19.6 and 27 GeV, while the z;2{FTPC} is larger 
than W2{4} in 39 GeV. This consistency can be also ob- 
served in Fig. H for 10 - 20% to 40 - 50% centrality bins. 
It indicates that the use of the first-order reaction plane 
(BBC event plane) to study the second harmonic flow 
eliminates flow fluctuations which are not correlated be- 
tween different harmonics. The first-order BBC reaction 
plane is struck by nucleon spectators for these beam en- 
ergies. The contribution of spectators makes the BBC 
event plane more sensitive to the reaction plane. This 
could partly explain the consistency between 1)2 {BBC} 
and i'2{4} mentioned above. More studies of the colli- 
sion energy dependence of non-flow and flow fluctuations 
will be discussed in another paper. 
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FIG. 7: (Color online) The «2 over e (CGC) as a function of pr for various collision centralities (10 - 20%, 30 - 40% and 
50-60%) in Au + Au collisions at midrapidity. Panel (a), (b), (c), (d) and (e) show the results for s/snn = 7.7, 11.5, 19.6, 
27 and 39 GeV respectively. The data are from U2{EtaSubs}. The error bars and shaded boxes represent the statistical and 
systematic uncertainties respectively, as described in Sec. IV C. 



C. Systematic uncertainties 

Different V2 methods show different sensitivities to 
non-flow correlations and V2 fluctuations. In previ- 
ous STAR publications, the differences between differ- 
ent methods were regarded as systematic uncertain- 
ties [nl, [3 ■ A great deal of progress has revealed that 
some of these differences are not due to systematic uncer- 
tainties in different methods, but due to different sensitiv- 
ities to non-flow and flow fluctuation effects [s^l . The 
four particle cumulant method is less sensitive to non- 
flow correlations [13, Ull and has a negative contribution 
from flow fluctuations. V2 measurements from the two 
particle cumulant method and the event plane method 
(the second harmonic event plane) have positive contri- 
butions from flow fluctuations as well as non-flow. It was 
also noticed that four particle cumulant results should 
be very close to flow in the reaction plane, while the 
two particle cumulant measures flow in the participant 
plane [H, Further, because of the large pseudora- 
pidity gap between the BBC/FTPC and TPC, 'y2{BBC} 
and i;2{FTPC} are most insensitive to non-flow correla- 
tions. 

We estimate the systematic uncertainty on event plane 
flattening methods for ?;2{EP} and z;2{EtaSubs} by the 
difference between them and find it to be negligible (be- 
low 1%). A 5% systematic uncertainty on W2{BBC}, 
W2{FTPC}, W2{EP} and ?;2{EtaSubs} is estimated by 
varying cut parameters (e.g. collision vertex position. 



the distance of closest approach to the primary vertex 
for the tracks, and the number of fit points used for re- 
construction of the tracks). The systematic uncertainties 
on V2{2} and W2{4} are based on the difference between 
Q-cumulant and GF-cumulant methods (5%) as well as 
cut variations (5%). All the percentage uncertainties are 
relative to the V2 value. 



V. DISCUSSION 

A. Transverse momentum and centrality 
dependence of V2 

The centrality dependence of pt differential V2 with 
respect to the initial eccentricity has been studied in de- 
tail for Au -I- Au and Cu -I- Cu collisions in ^/snn = 200 
and 62.4 GeV The larger magnitude of V2 in the 

more peripheral collisions could be due to the larger ini- 
tial eccentricity in coordinate space for the more periph- 
eral collisions. The participant eccentricity is the initial 
configuration space eccentricity of the participant nucle- 
ons defined by Eq. (3). The root- mean-square partici- 
pant eccentricity, Epart {2} As calculated from the Monte 
Carlo Glauber model jsa Isil (^Tab. HD and Color Glass 
Condensate (CGC) nrodel [llll (Tab.|lll|. The event 
plane is constructed from hadrons which have their ori- 
gin in participant nucleons. At the same time, the event 
plane resolution (i] sub-event) is less than 0.5. Thus, 
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FIG. 8: (Color online) Panel (a) shows the i)2{EP} vs. ri for 10 - 40% centrality in Au + Au collisions at i/sjvjv = 7.7, 11.5, 
19.6, 27, 39, 62.4 and 200 GeV. Panel (c) shows the ratio of V2 vs. rj for all ^/sjvjv with respect to the fit curve. Panel (b) shows 
the ii2{EP} vs. 'q/yhc!>.ra- Panel (d) shows the ratio of 112 vs. 7?/j/bcam for all ^smn with respect to the fit curve. The data for 
^ snn = 62.4 and 200 GeV are from refs. [l^. [61I. . The dashed red curves show the empirical fit to the results from Au + 
Au collisions at ^/snn = 7.7 GeV. The bands show the systematic uncertainties as described in Sec. IV C. 



what wc actually measure is the root-mean-square of V2 
with respeet to the partieipant plane [s^ . In this case, 
£part{2} is the appropriate measure of the initial geomet- 
ric anisotropy taking the event-by-event fluctuations into 
account [Hlf, [s^ [gO] ■ In Figs. |6] and [3 the centrality de- 
pendence of pt differential V2 over eccentricity is shown 
for Au + Au collisions at ^/s^vw" = 7.7, 11.5, 19.6, 27 
and 39 GeV. For all five collision energies, the centrality 
dependence of V2{pt) is observed to be similar to that at 
higher collision energies (62.4 and 200 GeV) of Au + Au 
and Cu -I- Cu colliding systems. That central collisions 
in general have higher V2/e than peripheral collisions is 
consistent with the picture that collective interactions arc 
stronger in collisions with larger numbers of participants. 



B. Pseudorapidity dependence 

The panel (a) of Fig. [5] shows V2 as a function of pseu- 
dorapidity for Au -f Au collisions at ^/Jnn = 7.7, 11.5, 
19.6, 27, 39, 62.4 and 200 GeV in mid-central (10 - 40%) 
collisions. The data for ^Jsmn = 62.4 and 200 GeV are 
from refs. [l^, |6ll, [s^] ■ To facilitate comparison with 62.4 
and 200 GeV data, the results of W2{EP} are selected for 



the rest of the collision energies. The 7.7 GeV data are 
empirically fit by the following function: 

V2{v) = Po + Pi'if + P2'r]'^ (45) 

with parameters po = 0.0450 ± 0.0002, pi = -0.0064 ± 
0.0015, p2 = -0.0024 ± 0.0017. For clarity, the panel (c) 
of Fig. [5] shows the ratio of W2(?7) with respect to this fit 
function. The pseudorapidity dependence of V2 indicates 
a change in shape as we move from ^/s^vw" = 200 GcV to 
7.7 GeV within our measured range -1 < 77 < 1. 

To investigate the collision energy dependence of the 
V2{r]) shape, in panel (b) and (d) of Fig. |S1 the same 
V2 results have been plotted as a function of pseudora- 
pidity divided by beam rapidity. The data of 7.7 GeV 
are fit by Eq. ^ with parameters po = 0.0450 ± 0.0002, 
Pl = -0.0279 ± 0.0064, p2 = -0.0464 ± 0.0325. The beam 
rapidities are 2.09, 2.50, 3.04, 3.36, 3.73, 4.20 and 5.36 
for yi]^ = 7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV 
respectively. After dividing pseudorapidity by the beam 
rapidity, the shape of V2 seems similar at all collision en- 
ergies. The approximate beam rapidity scaling on the 
V2{ri) shape suggests the change in shape may be related 
to the final particle density. Higher particle density indi- 
cates higher probability of interaction which can generate 
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FIG. 9: (Color online) The top panels show i'2{4} vs. Pt at midrapidity for various collision energies {^/smf = 7.7 GeV to 2.76 
TeV). The results for s/snn = 7.7 to 200 GeV are for Au + Au collisions and those for 2.76 TeV are for Pb + Pb collisions. 
The dashed red curves show the empirical fits to the results from Au + Au collisions at ^/snn = 200 GeV. The bottom panels 
show the ratio of U2{4} vs. pr for all ^/snn with respect to the fit curve. The results are shown for three collision centrality 
classes: 10 - 20% (al), 20 - 30% (bl) and 30 - 40% (cl). Error bars are shown only for the statistical uncertainties. 
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FIG. 10: (Color online) The experimental data (symbols) are the same as in Fig. |9] (b2). The lines represent the viscous 
hydrodynamic calculations from Ref. based on (a) MC-Glauber initial conditions and rj/s = 0.08 (b) MC-KLN initial 
conditions and rj/s = 0.20. 

larger collective flow. C. Energy dependence 

One of the most important experimental observations 
at RHIC is the significant V2 signal in the top energy of 
Au + Au collisions [l^ (more than 50% larger than at 
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FIG. 11: (Color online) The top panels show i;2{EtaSubs} vs. pt at midrapidity for various collision energies (^/sjviv" = 7.7 
GeV to 39 GeV). The bottom panels show the ratio of ii2{EtaSubs} vs. pr for all ^/snn with respect to the 39 GeV data. The 
results are shown for three collision centrality classes: 10-20% (al), 20-30% (bl) and 30-40% (cl). Error bars are shown 
only for the statistical uncertainties. 



the SPS [S^l)- It could be interpreted as the observation 
of a higher degree of thernialization than at lower colli- 
sion energies Q. The BES data from the RHIC-STAR 
experiment offers an opportunity to study the collision 
energy dependence of V2 using a wide acceptance detec- 
tor at midrapidity. Figure |9] shows the px dependence of 
V2{'i} from = 7.7 GeV to 2.76 TeV in 10-20% (al), 

20 - 30% (bl) and 30 - 40% (cl) centrality bins, where 
the ALICE results in Pb -I- Pb collisions at ^/s]v^'" = 
2.76 TeV are taken from Ref. [l^. The reasons to select 
the results of {4} for the comparison are the following: 
1) keep the method for V2 measurements consistent with 
the published results of ALICE; 2) W2{4} is insensitive to 
non-flow correlations. The 200 GeV data is empirically 
fit by a fifth order polynomial function. The parameters 
for the fit function are listed in Table IIVI For compari- 
son, the V2 from other energies are divided by the fit and 
shown in the lower panels of Fig. O We choose 200 GeV 
data as the reference because the statistical errors are 
smallest. For px below 2 GeV/c, the V2 values rise with 
increasing collision energy. Beyond = 2 GeV/c the V2 
results show comparable values within statistical errors. 

The increase of V2{pt) as a function of energy could be 
due to the ch ang e of chemical composition from low to 
high energies [22] and/or larger collectivity at the higher 
collision energy. The baryonic chemical potential varies a 
lot (20 - 400 MeV) from 200 to 7.7 GeV [H. The baryon 
over meson ratio is larger in lower collisions energies. The 
difference of V2 for baryon and meson, for example proton 



V2 < pion V2 for px below 2 GeV/c, could partly explain 
the collision energy dependence. Further, in Fig. [10] we 
compare the experimental data from Fig. ^ (b2) to the 
viscous hydrodynamic calculations As the collision 
energy varies from -^snn = 7.7 to 2760 GeV, the exper- 
imental data show larger splitting in the lower px region 
and converge at the intermediate range {px ~ 2 GeV/c); 
while, in the pure viscous hydrodynamic simulations, the 
splitting increases with px- The px dependence of the 
V2 ratio cannot be reproduced by pure viscous hydro- 
dynamic simulations with a constant shear viscosity to 
entropy density ratio (rj/s), and zero net baryon density. 
The comparison suggests that a quantitative study at 
lower collision energies requires a more serious theoret- 
ical approach, like 3D viscous hydro -I- UrQMD with a 
consistent equation of state at non-zero baryon chemical 
potential. 

Figure [TT] shows the energy dependence of 
?;2{EtaSubs}. Larger W2{EtaSubs} values are ob- 
served at higher collision energy for a selected px bin, 
but the Px dependence of the difference is quite different 
from W2{4}. The ratios to 39 GeV data for each collision 
energy first decrease as a function of px, then slightly 
increase in the px region of 1 - 2.5 GeV/c. The different 
trend of the energy dependence of V2 from 1^2 {4} and 
W2{EtaSubs} is interpreted as due to the different 
sensitivity of the V2 methods to non-fiow and/or fiow 
fiuctuations. 
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FIG. 12: (Color online) The U2{4} as a function of pr for 20 - 30% Au + Au collisions at = 7.7, 11.5, 39 and 200 GeV 

compared to corresponding results from UrQMD, AMPT default version, and AMPT with string melting version (3 and 10 
mb). The shaded boxes show the systematic uncertainties for the experimental data of 7.7, 11.5 and 39 GeV. The bottom 
panels show the ratio of data to the fit results of the models. 



D. Model comparisons 

To investigate the partonic and hadronic contribution 
to the final V2 results from different collision energies, 
transport model calculations from AMPT default (ver. 
1.11), AMPT string-melting (ver. 2.11) [H and UrQMD 
(ver 2.3) [65| are compared with the new data presented. 
The initial-parameter settings for the models follow the 
recommendation in the cited references. The AMPT de- 
fault and UrQMD models only take the hadronic interac- 
tions into consideration, while the AMPT string-melting 
version incorporates both partonic and hadronic interac- 
tions. The larger the parton cross section, the later the 
hadron cascade starts. 

Figure[T2lshows the comparison of pt differential f2{4} 
between model and data in the 20 - 30% centrality bin. 



The 200 GeV data are taken from [62| . The figure shows 
that UrQMD underpredicts the measurements at ^sjvw 
= 39 and 200 GeV in the pr range studied. The differ- 
ences are reduced as the collision energy decreases. That 
the ratio of data to UrQMD results are closer to 1 at the 
lower collision energy indicates that the contribution of 
hadronic interactions becomes more significant at lower 
collision energies. The AMPT model with default set- 
tings underpredicts the 200 GeV data, while the ratios of 
data to AMPT default results show no significant change 
from 7.7 to 39 GeV. The inconsistency between AMPT 
default and UrQMD makes the conclusion model depen- 
dent. The AMPT model with string- melting version with 
3 and 10 mb parton cross sections overpredicts the results 



at all collision energies from 7.7 to 200 GeV. A larger 
parton cross section means stronger partonic interactions 
which translate into a larger magnitude of V2 ■ The differ- 
ence between data and these AMPT model calculations 
seems to show no significantly systematic change vs. col- 
lision energies. However, a recent study with the AMPT 
model suggests hadronic potentials affect the final V2 re- 
sults significantly when the collision energy is less than 
^/sWn = 39 GeV H^. 



VI. SUMMARY 

We have presented elliptic flow, V2, measurements from 
Au + Au collisions at ^/s/viv = 7.7, 11.5, 19.6, 27 and 
39 GeV for inclusive charged hadrons at midrapidity. To 
investigate non-flow correlations and V2 fluctuations, var- 
ious measurement methods have been used in the anal- 
ysis. The difference between i'2{2} and V2{4:} decreases 
with decrease in collision energy, indicating that non-flow 
contribution and/or flow fluctuations decrease with a de- 
crease in collision energy. The centrality and pT depen- 
dence of V2 are similar to that observed at higher RHIC 
collision energies. A larger V2 is observed in more pe- 
ripheral collisions. The pseudorapidity dependence of V2 
indicates a change in shape from 200 GeV to 7.7 GeV 
within the measured range -1 < ij < 1, but the results of 
V2 versus pseudorapidity scaled by beam rapidity shows 
a similar trend for all collision energies. The compari- 
son with Au -I- Au collisions at higher energies at RHIC 
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(yi^ = 62.4 and 200 GcV) and at LHC (Pb + Pb colli- 
sions at ^ysffN = 2.76 TcV) shows the '(;2{4} values at low 
Pt {pt < 2.0 GcV/c) increase with increase in collision 
energy implying an increase of collectivity. The current 
viscous hydrodynamic simulations cannot reproduce the 
trend of the energy dependence of V2{pt)- 

The agreement between the data and UrQMD, which 
is based on hadronic rescatterings, improves at lower 
collision energies, consistent with an increasing role of 
the hadronic stage at these energies. The inconsistency 
between AMPT default and UrQMD makes the conclu- 
sion model dependent. The comparison to AMPT model 
calculations seems to show no significantly systematic 
change vs. collision energy, but improved calculations 
including harmonic potentials may change the V2 values 
from AMPT models at lower collision energies. 

These results set the baseline to study the number of 
constituent quark scaling of identified hadron V2 ■ It also 
sets the stage for understanding the collision energy de- 
pendence of V2 in the regime where the relative contribu- 



tion of baryon and mesons vary significantly. 
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TABLE II: Summary of centrality bins, average number of participants (A'^part), number of binary collisions (iVcoii), reaction 
plane eccentricity (erp), participant eccentricity (epart), root-mean-square the participant eccentricity epart{2} and transverse 
area (Spart) from MC Glauber simulations at ^/snn ~ 7.7, 11.5, 19.6, 27, and 39 GeV. The errors are systematic uncertainties. 
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0.32 ± 0.04 


0.37 ± 0.03 


0.39 ± 0.03 


12.6 ± 0.8 


/10 PiOPT, 


7Q -L 1 
( ± lU 


98 ± 18 


0.36 ± 0.04 


0.43 ± 0.03 


0.46 ± 0.03 


10.1 ± 0.9 


OU-DU /O 


AA -L Q 


52 ± 14 


0.39 ± 0.04 


0.50 ± 0.03 


0.53 ± 0.03 


7.8 ± 1.0 


fio 70*?;^^ 

DU- / U /O 


9(^-1-7 
ZD ± / 


25 ± 9 


0.40 ± 0.05 


0.58 ± 0.04 


0.62 ± 0.04 


5.8 ± 1.1 


70 80*?;^^ 

/ U-oU /O 


1/1-1-1^ 
1^ ± D 


12 ± 6 


0.37 ± 0.06 


0.68 ± 0.05 


0.71 ± 0.05 


3.7 ± 1.2 








Au 4- Au at ,ysjvjv - 


= 19.6 GeV 






U-070 


338 ± 2 


800 ± 27 


0.044 ± 0.006 


0.102 ± 0.003 


0.117 ± 0.003 


25.6 ± 0.4 


0- 1 U /o 


OQO _i_ 

zoy ± D 


643 ± 20 


0.11 ± 0.01 


0.15 ± 0.01 


0.16 ± 0.01 


23.0 ± 0.3 


1 90% 


99c; -L Q 

zzo ± y 


458 ± 24 


0.18 ± 0.02 


0.22 ± 0.02 


0.24 ± 0.02 


19.5 ± 0.5 


90 '^0% 


1 p;e 1 in 
100 ± lU 


284 ± 26 


0.26 ± 0.03 


0.30 ± 0.02 


0.32 ± 0.02 


15.6 ± 0.7 


Qfi A0% 


1 08 -1- 10 

lUo it lU 


170 ± 23 


0.32 ± 0.04 


0.37 ± 0.03 


0.40 ± 0.03 


12.5 ± 0.8 


40 ^0*5^ 


71 -1- 1 
/ 1 it lU 


96 ± 18 


0.36 ± 0.04 


0.43 ± 0.03 


0.46 ± 0.03 


10.0 ± 0.9 


^0 RO*J^ 

OU-DU /O 


44 4- Q 
H:'^ it y 


51 ± 13 


0.39 ± 0.04 


0.50 ± 0.03 


0.53 ± 0.03 


7.8 ± 1.0 


flO 70% 

DU- i U /O 


9^-1-7 

ZO it ( 


25 ± 8 


0.40 ± 0.05 


0.58 ± 0.04 


0.62 ± 0.04 


5.8 ± 1.1 


70 S0% 
/ u-ou /o 


14 -1- 
1^ it 


12 ± 5 


0.37 ± 0.06 


0.68 ± 0.05 


0.71 ± 0.05 


3.7 ± 1.2 








Au -f Au at \/smn 


= 27 GeV 






U-070 


343 ± 2 


841 ± 28 


0.040 ± 0.005 


0.100 ± 0.002 


0.114 ± 0.003 


25.8 ± 0.4 


o-lU/o 


zyy lb D 


694 ± 22 


0.10 ± 0.01 


0.14 ± 0.01 


0.16 ± 0.01 


23.4 ± 0.3 


1 90% 
lu-zu /o 


9QQ -L Q 

ZOO ± y 


497 ± 26 


0.18 ± 0.02 


0.21 ± 0.02 


0.23 ± 0.02 


19.8 ± 0.5 


90 '^0% 
ZU-OU /o 


IDD it 11 


312 ± 28 


0.26 ± 0.03 


0.29 ± 0.02 


0.32 ± 0.02 


15.9 ± 0.7 


•^0 /10% 

ou—^u /o 


11/1 -U 11 
1 i-y: it 11 


188 ± 25 


0.32 ± 0.04 


0.37 ± 0.03 


0.39 ± 0.03 


12.8 ± 0.9 


'iU-OU /o 


7p; -L 1 
/ ± lU 


106 ± 20 


0.37 ± 0.04 


0.43 ± 0.03 


0.46 ± 0.03 


10.2 ± 0.9 


PiO R0% 
OU-DU /O 


A7 -L Q 


56 ± 15 


0.39 ± 0.05 


0.50 ± 0.03 


0.53 ± 0.03 


7.9 ± 1.0 


fiO 70% 
DU- i U /O 


97 -1- 8 
Z / ± 


27 ± 10 


0.40 ± 0.05 


0.58 ± 0.05 


0.61 ± 0.05 


5.8 ± 1.2 


70 80% 
/ U-oU /O 


1 /I -1- 

1^ it D 


12 ± 6 


0.37 ± 0.06 


0.68 ± 0.05 


0.71 ± 0.05 


3.6 ± 1.3 








Au + Au at ,ysMN 


= 39 GeV 






0-5% 


342 ± 2 


853 ± 27 


0.042 ± 0.006 


0.101 ± 0.003 


0.115 ± 0.003 


25.9 ± 0.4 


5-10% 


294 ± 6 


687 ± 21 


0.10 ± 0.01 


0.14 ± 0.01 


0.16 ± 0.01 


23.3 ± 0.3 


10-20% 


230 ± 9 


492 ± 26 


0.18 ± 0.02 


0.21 ± 0.02 


0.23 ± 0.02 


19.8 ± 0.5 


20-30% 


162 ± 10 


306 ± 27 


0.26 ± 0.03 


0.30 ± 0.02 


0.32 ± 0.02 


16.0 ± 0.7 


30-40% 


111 ± 11 


183 ± 24 


0.32 ± 0.04 


0.37 ± 0.03 


0.39 ± 0.03 


12.8 ± 0.8 


40-50% 


74 ± 10 


104 ± 20 


0.36 ± 0.04 


0.43 ± 0.03 


0.46 ± 0.03 


10.3 ± 1.0 


50-60% 


46 ± 9 


55 ± 14 


0.39 ± 0.04 


0.50 ± 0.03 


0.53 ± 0.03 


8.0 ± 1.0 


60-70% 


26 ± 7 


27 ± 9 


0.40 ± 0.04 


0.58 ± 0.04 


0.61 ± 0.04 


5.9 ± 1.1 


70-80% 


14 ± 5 


12 ± 6 


0.37 ± 0.05 


0.67 ± 0.05 


0.71 ± 0.05 


3.8 ± 1.2 
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TABLE III: The epart{2} and transverse area (Spart) from the Color Glass Condensate (CGC) model [ssl - ls^ calculations in 
Au + An collisions at ^snn ~ 7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV. The errors are systematic uncertainties. 



Centrality (%) 


0-5% 


5 - 10% 


10 - 20% 20 - 30% 


30 - 40% 40 - 50% 


50 - 60% 


60 - 70% 70 - 80% 








All 1 An Q t / Q HT-r^r 


— 7 7 HeV 






epart{2} 
\ Impart/ ^illi j 


0.104 ±0.005 
95 Q + 1 3 


0.19 ±0.01 
21.8 ±1.1 


0.29 ±0.01 0.39 ±0.02 
175 + nq 134 + 07 


0.47 ±0.02 0.54 ±0.03 
109 + 05 77 + 04 


0.59 ±0.03 
5.5 ± 0.3 


0.62 ±0.03 0.51 ±0.02 
36+09 18+01 








All —I— All / Q HTM 


— 11 5 GeV 






£part{2} 
\ Impart/ ) 


0.104 ±0.005 
959+1 9 


0.19 ±0.01 
21.2 ±1.1 


0.29 ±0.01 0.39 ±0.02 
i7n + nQ i3n + n7 


0.47 ±0.02 0.53 ±0.03 
QQ+05 75+04 

C7.C/ III KJ'O 1 'O III 


0.59 ±0.03 
5.4 ± 0.3 


0.62 + 0.03 0.51 ±0.02 
35+09 18+01 








All -I— All / Q ATM 


— 1 Q 6 CpV 






epart{2} 
\ Impart/ J 


0.105 ±0.005 
24.4 ± 1.2 


0.19 ±0.01 
90 fi + 1 n 


0.29 ±0.01 0.39 ±0.02 


0.47 ±0.02 0.53 ±0.03 
Q7+05 73+04 

iJ . 1 III KJ.O 1 .O III 


0.58 ±0.03 
5.3 ± 0.3 


0.61 ±0.03 0.51 ±0.02 
35+09 18+01 








All 1 All £a4" / c AT A 


- = 27 GeV 






epart{2} 


0.105 ±0.005 
24.1 ±1.2 


0.19 ±0.01 

90 3 + 1 n 


0.29 ±0.01 0.39 ±0.02 

164 + nS 195 + 06 


0.47 ±0.02 0.53 ±0.03 

Qfi+05 79+04 


0.58 ±0.03 
5.3 ± 0.3 


0.61 + 0.03 0.51 ±0.02 
35+09 18+01 








All —1— All / Q AT A 


- — 39 QeV 






£part{2} 
(Spart) (fm^) 


0.105 ±0.005 
23.9 ± 1.2 


0.19 ±0.01 
20.1 ± 1.0 


0.29 ±0.01 0.39 ±0.02 
16.2 ±0.8 12.4 ±0.6 


0.47 ±0.02 0.53 ±0.03 
9.5 ±0.5 7.2 ±0.4 


0.58 ±0.03 
5.3 ±0.3 


0.61 + 0.03 0.50 ±0.02 
3.5 ±0.2 1.8 ±0.1 








Au + Au at ^/sNN 


= 62.4 GeV 






epart{2} 

(Spart) (fm^^) 


0.105 ±0.005 
23.7± 1.2 


0.19 ±0.01 
20.0 ±1.0 


0.29 ±0.01 0.39 ±0.02 
16.1 ±0.8 12.3 ±0.6 


0.47 ±0.02 0.53 ±0.03 
9.4 ±0.5 7.2 ±0.4 


0.58 ±0.03 
5.3 ±0.3 


0.61 ±0.03 0.50 ±0.02 
3.5 ±0.2 1.8 ±0.1 








Au + Au at ^ysNN 


= 200 GeV 






£part{2} 
(Spart) (W) 


0.104 ±0.005 
23.7 ± 1.2 


0.19 ±0.01 
20.0 ±1.0 


0.29 ±0.01 0.39 ±0.02 
16.1 ±0.8 12.3 ±0.6 


0.47 ±0.02 0.53 ±0.03 
9.4 ±0.5 7.2 ±0.4 


0.57 ±0.03 
5.3 ±0.3 


0.60 ±0.03 0.49 ±0.02 
3.6 ±0.2 1.9 ±0.1 



TABLE IV: Summary of the parameters for the fit functions to the results of 112 {4} vs. pr in Au+Au coUisions at ^/snn ~ 200 
GeV. 



Parameters 


Pa 


Pi 


P2 


pa 


P4 


P5 


10 - 20% 


-0.00730 ±0.00114 


0.10785 ±0.00598 


-0.03941 ±0.01038 


0.01508 ±0.00767 


-0.00411 ± 0.00246 


0.00041 ± 0.00028 


20 - 30% 


-0.00890 ± 0.00096 


0.14250 ±0.00500 


-0.05206 ± 0.00869 


0.02156 ± 0.00642 


-0.00685 ± 0.00206 


0.00077 ± 0.00023 


30 - 40% 


-0.00581 ±0.00206 


0.14526 ±0.01089 


-0.00529 ±0.01910 


-0.02409 ±0.01419 


0.00797 ± 0.00456 


-0.00084 ± 0.00052 



